Ab initio study of proper topological ferroelectricity in layered perovskite La2Ti207 
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We present a first-principles investigation of ferroelectricity in layered perovskite oxide La2Ti2 07 
(LTO), one of the compounds with highest Curie temperature known (1770 K). Our calculations 
reveal that LTO's ferroelectric transition results from the condensation of two soft modes that 
have the same symmetry and are strongly coupled anharmonically. Further, the leading instability 
mode essentially consists of rotations of the oxygen octahedra that are the basic building block 
of the perovskite structure; remarkably, because of the particular topology of the lattice, such Oe 
rotations give raise to a spontaneous polarization in LTO. The effects discussed thus constitute an 
example of how nano-structuring - provided here by the natural layering of LTO - makes it possible 
to obtain a significant polar character in structural distortions that are typically non-polar. We 
discuss the implications of our findings as regards the design of novel multifunctional materials. 
Indeed, the observed proper ferroelectricity driven by Oe rotations provides the ideal conditions to 
obtain strong magnetoelectric effects. 

PACS numbers; 77.84.-s, 61.50.Ah, 75.85. +t, 71.15.Mb 
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I. INTRODUCTION 

Because of their physical appeal and technological im- 
portance, ferroelectrics and related materials have been 
the object of continued attention for decades^i""^ Bulk 
oxides with the ideal perovskite structure ~ ranging 
from classic ferroelectric BaTiOs to strong dielectric 
Bai_a:Sra;Ti03, or from piezoelectric PbZri_2:Tia;03 to 
relaxor (PbMgi/3Nb2/303)i-:.-(PbTi03):. - have been 
especially well studied. Indeed, partly thanks to a num- 
ber of key contributions from first-principles theory,^- we 
now understand the fundamental atomistic origin of the 
ferroelectric (FE) and response properties of the most 
important members of this family. During the past 
decade, the focus has increasingly shifted towards nano- 
structured materials, especially in the form of thin films. 
Work on films has led to a better understanding of how 
elastic (i.e., the epitaxial strain exerted by a substrate) 
and electric (e.g., the imperfect screening of the depo- 
larizing field associated to particular metallic electrodes) 
boundary conditions affect the FE statci^ Further, it has 
been shown that exotic FE properties can be obtained in 
artificially created super-lattices, as e.g. in the recently 
discussed PbTi03/SrTi03 heterostructuresi^ 

The emergence of magnetoelectric multiferroics^ (i-C-, 
materials with coupled magnetic and FE orders) has 
contributed to refuel interest in ferroelectrics, especially 
in what regards unconventional mechanisms for ferro- 
electricity. If we restrict ourselves to oxides with the 
ideal ^4503 perovskite structure, ferroelectricity (usu- 
ally driven by a i3-site transition metal with an ndP elec- 
tronic configuration) and magnetism (which requires lo- 
calized d 01 f electrons) seem to be mutually exclusive, 
the known exceptions being very scarce.-S*2- So far, the 
most notable ways around this problem consist in (i) hav- 
ing ferroelectricity driven by the ^-site cation, as in room 
temperature multiferroic BiFe03,^'' and (ii) moving away 
from the ideal perovskite structure and resorting to other 
mechanisms for ferroelectricity, as e.g. in the improper 
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FIG. 1; (Color on-line.) Different views of the layered per- 
ovskite structure of the A„BnOsn+2 compounds for n = 4, 
where n is the number of perovskite-like planes within a layer. 
Only the A atoms (as balls) and BOg octahedra are shown. 
The defined Cartesian axes x, y, and z (which follow the con- 
vention adopted in Ref. 13^ ) are used throughout the paper. 



ferroelectrics YMnO^iiii^ and Ca3Mn207J^ 

This is the context of our work on the layered 
perovskite oxide La2Ti207 (LTO), whose structure is 
sketched in Fig. 1. LTO is one of the highest-temperature 
ferroelectrics known, with a Curie point (Tq) of approx- 
imately 1770 and is lately being considered for ap- 
plications as a high-T piezoelectric i^^'^^ The FE transi- 
tion temperature of LTO is enormous, especially when 
compared with the Tc's of cubic perovskite titanates one 
might have expected to be similar to LTO: most signif- 
icantly, BaTiOs (BTO) becomes ferroelectric at about 
400 K, and the related compound SrTi03 remains (a 
quantum) paraelectric down to K. In order to explain 
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the surprising behavior of LTO, we need to identify the 
atomistic origin of its FE instabihty. LTO might owe 
its very high Tq to the kind of mechanisms known to 
operate in the cubic titanates (in essence, long-range 
dipole-dipole interactions that destabihze the non-polar 
phasei^iiSi) , which might somehow be enhanced by the 
peculiar topology of the LTO lattice. Alternatively, LTO 
might present some new form of very strong FE instabil- 
ity. Either way, the study of LTO will provide us with 
information that might be useful for the design of new 
ferroelectrics. 

Let us also note that LTO is a member of the family of 
oxides with general formula La„Ti„03„4.2,^^ whose struc- 
ture can be seen as the ideal cubic perovskite periodically 
truncated along the [011]c direction of the cubic lattice, 
n being the number of perovskite-like planes within one 
layer (see Fig. 1). These structures are thus related to the 
well-known Ruddlesden-Popper, Aurivillius, and Dion- 
Jacobson families of layered perovskites,- ° for which the 
truncation direction is [001] c and which include very fa- 
mous members such as (La,Ba)2Cu04, the parent com- 
pound of high-Tc superconductors. The basic features of 
the electronic structure of the La„Ti„03„+2 compounds 
are controlled by n: The La and O atoms can be assumed 
to be in their most common ionization states - i.e., La'^^ 
and ~, which implies that the Ti cations will present 
a positive charge of 3 -I- 4/n. Accordingly, the number of 
3d electrons in the Ti atoms will be 1 — 4/n, which allows 
us to move quasi-continuously from the Ti-Sd" configu- 
ration of the n=4 compound (that is our La2Ti207, the 
family member with smallest n reported in the litera- 
ture) to the Ti-3o?^ configuration of the n—oo compound 
(which has the prototype perovskite structure). The va- 
riety of electronic phenomena observed for intermediate 
values of n - including phases of semi-conducting, normal 
metallic, and low-dimensional metallic character— - con- 
stitutes an additional motivation to study in detail the 
structural behavior of the relatively simple end member 
La2Ti207. 

The paper is organized as follows. In Section II we de- 
scribe the theoretical approach and first-principles meth- 
ods used in this work. We present and discuss our re- 
sults in Section III, which is split in the following way. 
In Section III. A we describe our results for the struc- 
ture of the high-temperature paraelectric (PE) and FE 
phases of LTO. In Section III.B we show that the PE 
phase presents a strong FE instability whose topologi- 
cal nature is discussed in detail. In Section III.C we 
describe the phase transition between the PE and FE 
phases, and discuss the energetics of the transformation. 
In Section III.D we present our results for the sponta- 
neous polarization, dielectric, and piezoelectric proper- 
ties. Having established LTO's peculiarity, which is made 
manifest by means of a detailed comparison with proto- 
type compound BTO, in Section III.E we show that LTO 
also presents some features that are clearly reminiscent 
of the usual FE oxides with the ideal perovskite struc- 
ture. In Section III.F we outline the exciting implications 



that our results have for the design of novel magneto- 
electric materials. Finally, in Section IV we summarize 
and present our conclusions. Whenever it is possible, 
we compare our first-principles results with the (scarce) 
experimental information available for LTO. 

II. METHODOLOGY 

A. Theoretical approach to La2Ti207 

We adopted the usual first-principles approach to the 
investigation of structural phase transitions of the dis- 
placive type, which is routinely applied with great success 
to FE perovskite oxidesi^ Here we describe a simplified 
version of such an approach that is sufficient for our study 
of LTO (i.e., we will assume one-dimensional instability 
modes, will not discuss how to deal with cell strains, etc.). 
In essence, one needs to identify a reference equilibrium 
phase of high symmetry (HS) - which corresponds to 
the high-temperature PE phase of the compound - and 
study its stability against all possible structural distor- 
tions. More precisely, one writes the energy of the crystal 
as the following Taylor series: 

^ = + IY1 Kn^nU^Un + 0{U^) , (1) 

rri.n 

where is the energy of the HS phase. The Um vari- 
ables represent the structural distortions of the refer- 
ence configuration; for the study of a FE transition, one 
can typically restrict oneself to the distortions compati- 
ble with the reference unit cell, i.e., those associated to 
the F-point of the Brillouin zone of the PE phase. The 
so-called force-constant matrix K is the central quantity 
one needs to compute, as its negative eigenvalues (if any) 
correspond to unstable structural distortions - i.e., soft 
modes - that may result in a phase transition. Let us 
use the term mode stiffness to refer to the eigenvalues 
of which we denote by with s running from 1 to 
the dimension of the force-constant matrix. Note that 
the magnitude of a negative Ks determines the strength 
of the structural instability, and thus the likelihood of 
observing it experimentally. (In general, one may find 
several, possibly competing, instabilities of a HS phase; 
it is by no means guaranteed that all of them will lead 
to experimentally observable phase transitions.) Once a 
soft mode is identified, one can readily study the corre- 
sponding low-symmetry (LS) phase - by distorting the 
crystal according to the soft mode eigenvector and re- 
laxing the resulting structure - and the energetics of the 
instability - by computing the usual double-well poten- 
tial connecting the HS and LS phases. 

In order to apply this program to LTO, we had to 
tackle one fundamental difficulty: We could not find any 
experimental information on the structure of the high- 
temperature PE phase of this compound.-- Note that 
this is never a problem when one works with the usual FE 
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FIG. 2: (Color on-line.) Structures of the Cmc2i (FE) and 
Cmcm (PE) phases of La2Ti2 07 studied in this work. The 
small (red), medium (blue), and big (violet) balls represent 
O, Ti, and La atoms, respectively. The shaded polyhedra are 
top- viewed Oe groups. The two phases have the same 44-atom 
conventional cell depicted in the figure; the 22-atom unit cell, 
which is also common to both, is indicted with dashed lines 
in the Cmcm case. The defined Cartesian axes x, y, and 
z (which follow the convention adopted in Ref. [3^ are used 
throughout the paper. The symmetry-independent atoms are 
labeled as in Table I; note that oxygens 0(1) and 0(3) of 
the Cmcm structure split, respectively, into 0(l)/0(2) and 
0(3)/0(4) in the Cmc2i structure. 



perovskites, where the ideal cubic perovskite structure is 
the reference phase of choice. Such a choice is obviously 
correct in cases like those of BTO or PbTiOs, where the 
cubic PE phase is experimentally accessible for T > Tc; 
more remarkably, this choice is also the most physically 
sound one for materials (e.g., BiFeOs) whose cubic PE 
phase is not easy to access experimentally, as the samples 
tend to melt before reaching the corresponding transition 
temperature. In the general case, the problem of choos- 
ing an appropriate HS reference phase for the theoret- 
ical study of displacive phase transitions has long been 
solved. The basic idea is to look for pseudo-symmetries 
(i.e., slightly broken symmetries) of the known LS struc- 
ture; we can thus identify possible HS phases that would 
transform into the known LS structure upon a relatively 
small distortion. This is a very powerful strategy that 
can lead to the discovery of complex phase transition se- 
quences (e.g., when more than one possible HS phases 
are found) and, in particular, has been used to identify 
previously unnoticed FE transitions j^^^ 

Using widely available crystallographic tools,— we ap- 
plied the pseudo-symmetry analysis to LTO and obtained 
a very clear prediction: We found that the FE phase of 
this compound, which presents space group C'mc2i and 



a 22-atom primitive unit cell, is most likely associated to 
a PE phase with the same primitive cell and space group 
Cmcm (see sketch in Fig. 2). We thus performed our 
first- principles study using this Cmcm phase as our HS 
reference structure. 

A few additional points are in order: (1) The pseudo- 
symmetry analysis resulted in relatively large atomic dis- 
placements connecting the HS and LS phases, with max- 
imum values of about 0.4 A corresponding to the La 
atoms. Note that the magnitude of such distortions 
is expected to reflect the associated HS-LS transition 
temperaturCf^^ which is indeed very high in this case. (2) 
Experimental studies of the compounds Sr2Nb207 and 
Ca2Nb207^^ which share with LTO the same layered 
structure and a similar (ndP) electronic configuration of 
the transition metal atoms, show that they present a 
high-temperature phase with Cmcm space group. (3) 
The choice of Cmcm as our PE space group results in spe- 
cific predictions for the symmetry of the soft modes that 
would be associated with a FE transition to the Cmc2i 
phase. Indeed, the leading instability should transform 
with the Biu irreducible representation of mmm, the 
point group of the Cmcm phase. As we will see, our 
first- principles results confirmed these expectations, thus 
supporting our choice of space group for the PE phase. 



B. Details of the calculations 

We used the local density approximation (LDA) to 
density functional theory (DFT) as implemented in the 
first-principles package VASPi^ To represent the ionic 
cores we used the projector-augmented wave schemCf^ 
solving explicitly for the following electrons: La's 5p, 5d, 
and 6s; Ti's 3p, 3d, and 4s; and O's 2s and 2p. The elec- 
tronic wave functions were represented in a plane-wave 
basis truncated at 400 eV. We always worked with the 44- 
atom cell of LTO sketched in Fig. 2 (which is the conven- 
tional cell for both Cmcm and Cmc2i phases), and used 
a 6x 1 x5 fc-point grid for Brillouin zone integrations. We 
checked these calculation conditions were well-converged 
by monitoring the computed equilibrium structure, bulk 
modulus, and phonon frequencies of the Cmcm phase. 
For the calculation of the force-constant matrix K and 
the dielectric and piezoelectric responses, we employed a 
simple finite-displacement scheme and the corresponding 
linear-response tensor formulas, which can be found e.g. 
in Ref. l31j. We only computed the lattice-mediated part 
of the static dielectric response, which has been repeat- 
edly shown to dominate the effect in ferroelectric oxides. 

Our BTO simulations were analogous to the above de- 
scribed ones. We solved explicitly for Ba's 5s, 5p, and 
6s electrons (treating Ti and O as above), and used a 
400 eV cut-off for the plane-wave basis. We worked with 
the 5-atom unit cell of BTO, and used a 6x6x6 /c-point 
grid for Brillouin zone integrations. 
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III. RESULTS AND DISCUSSION 

A. Structure of the Cmcrn and Cmc2i phases 

Table I shows our results for the equilibrium struc- 
ture of the Cmcrn (PE) and Cmc2i (FE) phases of LTO 
considered in this study. For Cmc2i we also report the 
experimental structure measured at 1173 K by Ishizawa 
et al— using X-ray diffraction. 

While the overall agreement between the experimen- 
tal and theoretical FE structures is acceptable, the de- 
viations affecting some atomic positions and lattice con- 
stants are larger than what is usual for this type of calcu- 
lations. For example, the predicted position of the La(l) 
atom differs notably from the experimentally determined 
one; more significantly, the computed a and c lattice con- 
stants deviate from the experimental values by almost a 
3%. We attribute these discrepancies to the fact that 
we are comparing our first-principles results, which cor- 
respond to the limit of K, with structural data taken at 
very high temperatures. Given that thermal expansion 
and other temperature-driven effects are not included in 
our simulations, our results seem compatible with the 
experimental information. 

B. Ferroelectric instabilities of the Cmcrn phase 

As described in Section II. A, we studied the structural 
stability of the Cmcm phase against F-point distortions 
(compatible with the PE unit cell) by computing the cor- 
responding force-constant matrix K. From the diagonal- 
ization of K we obtained two negative eigenvalues that 
correspond to two structural instabilities. The computed 
mode stiffnesses are —0.81 eV/A^ and —0.01 eV/A^, re- 
spectively. Hence, according to our calculations, LTO's 
Cmcm phase presents a structural instability comparable 
in strength with the FE soft mode of BTO's PE phase (for 
which we obtained —2.74 eV/A^), as well as a marginally 
unstable mode with nearly zero stiffness. 

Both instabilities transform with the _Bi„ irreducible 
representation of the mmm point group of the PE phase; 
more precisely, they are infra-red active (polar) modes 
that involve the development of a polarization along the 
z direction defined in Fig. 2. Such _Bi„ modes break the 
C21/, C2xi I, and az point symmetries of the PE phase 
(where / is the spatial inversion, C2a stands for a two- 
fold rotation around axis a, and ctq is a mirror plane per- 
pendicular to direction a), leading to the Cmc2i space 
group. Hence, the obtained Biu soft modes are exactly 
the kind of instabilities that can drive a ferroelectric 
phase transition between the Cmcm and Cmc2i phases. 
Our first-principles results thus support the correctness 
of our working hypothesis, i.e., that the studied Cmcm 
structure is indeed the high-symmetry PE phase of LTO. 

We can gain insight into the origin of ferroelectricity in 
LTO by inspecting the eigenvector of the strongest insta- 
bility mode, denoted by in the following. In essence. 



TABLE I: Computed equilibrium structures of the Cmcm and 
Cmc2i phases of La2Ti207 discussed in the text. We show 
in parenthesis the experimental values reported in Ref. for 
the Cmc2i phase. 
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involves a rotation of the Og oxygen octahedra around 
the X axis, as sketched in Fig. 3(a). Hence, the displace- 
ments of the equatorial oxygens amount to most of the 
eigenvector (11% of the norm of is associated to 0(1) 
displacements, 51% to 0(3), and 28% to 0(5), using the 
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FIG. 3: (Color on-line.) Panel (a): sketch of the largest 
atomic displacements associated to the strongest instability 
mode (^j^) obtained for the Cmcm phase of La2Ti207. We 
show the displacements corresponding to one layer composed 
of n = 4 perovskite-like planes. The arrows on the side rep- 
resent the electric dipoles associated to the displacement of 
oxygens in different y-planes (see text). The dipoles are la- 
beled as the oxygen atoms in Fig. 2. Panel (b): sketch of a 
typical anti-ferrodistortive mode occurring in an ideal (non- 
layered) perovskite structure. 



labels defined in Fig. 2 for the oxygen atoms); there are 
also significant La displacements (9% of the norm) , while 
the Ti atoms and apical oxygens 0(6) and 0(7) have a 
negligible participation in The character of the sec- 
ond soft mode (^2) niore complex: It is dominated by 
the displacements of oxygens 0(1) and 0(5) (with 21% 
and 35% of the norm, respectively) and involves a signif- 
icant deformation of the Og octahedra; it also presents 
a large participation of the La atoms (about 24% of the 
norm). 

These FE instabilities are very different from the ones 
that are usual among ABO^ oxides with the ideal per- 
ovskite structure. The inset of Fig. 4(b) shows the rep- 
resentative case of BTO: the Ti cation moves away from 
the center of the Og octahedron, which is only slightly 
distorted, giving raise to a large electric dipole. Fur- 
ther, the FE instability in BTO and related materials is 
known to originate from the strong interactions between 
such local dipoles.^® In contrast, the dominant FE insta- 
bility found in LTO consists of Og octahedra rotations, 
the off-centering of the Ti atoms being negligible. As 
quantified in Section III.D, such a pattern of atomic dis- 
placements does not lead to a large local dipole, which 
suggests that the mechanisms responsible for the ferro- 
electricity in LTO have to be of a different nature. 

Interestingly, structural instabilities involving Og oc- 
tahedra rotations are very common among ABO3 per- 
ovskites, and are usually termed anti-ferrodistortive 
(AFD). Indeed, AFD modes as the one sketched in 
Fig. 3(b) are the driving force for most of the struc- 
tural phase transitions occurring in these compounds, 
the examples including crystals as well-known as SrTiOa, 
LriMnOs and LnNiOs (where Ln is a lanthanide), mul- 
tiferroics BiMOs (where M is a 3d transition metal), 
etc. Accordingly, there is extensive literature devoted 
to the study and classification AFD modes in ABO3 
perovskiteSf^^^ and it is known that size (e.g., the in- 
compatibility of the ionic radii of the A and B cations 



to form a cubic perovskite lattice) and chemical (as 
in the Bi-based compounds that display the so-called 
stereochemical activity) effects are usually responsible for 
the occurrence of AFD distortions f2S, Hence, it is not a 
surprise to find that layered perovskite structures may 
present AFD-like instabilities that allow for a better ful- 
fillment of steric and/or chemical constraints at a local 
level. It is not our goal here to discuss the occurrence and 
origin of such soft modes in LTO-like layered perovskites; 
such a study should probably involve consideration of a 
number of representative crystals (e.g., a few members 
of the AnBnO^n+2 family) and falls beyond the scope of 
this work. Suffice it to say that the leading FE instabil- 
ity that we found in LTO is essentially analogous to the 
AFD distortions that are ubiquitous among perovskite 
oxides. 

However, we must note a critical difference between 
LTO's AFD-like mode depicted in Fig. 3(a) and the AFD 
modes occurring in ideal perovskite structures [Fig. 3(b)]: 
the former causes a spontaneous polarization, while the 
latter do not. To better explain this, in Fig. 3(a) we in- 
dicate with arrows the electric dipoles that appear as a 
consequence of the displacement of the oxygen atoms fol- 
lowing the eigenvector. The oxygens in each y-oriented 
plane give raise to a dipole along the z direction. Using 
the definitions in Fig. 3(a), the total dipole associated to 
one layer would be p'^^'' = p^-\- 2p^ + 2pi . Since there is 
no symmetry relationship between the displacements of 
the 0(1), 0(3), and 0(5) oxygens in the eigenvector, it 
follows immediately that p^^y^'^ will be different from zero. 
Moreover, even if oxygens of different types were to dis- 
place by the same amount, so that p^ = — pg = Pi = p, 
we would still have p'^y®'' = p ^ 0, as the number of 
oxygen planes in each layer is odd. Finally, note that 
LTO's unstable mode involves an identical distortion 
of all layers in the structure, and thus gives raise to a net 
macroscopic polarization. 

On the other hand, AFD distortions of the ideal (non- 
layered) perovskite structure do not result in a macro- 
scopic polarization. In Fig. 3(b) we show a pattern of 
Og rotations around the indicated ; the situation 

is very similar to the one depicted in Fig. 3(a), except 
that in this case the network of Og octahedra is not 
truncated. Again, we indicate with arrows the electric 
dipoles that appear as a consequence of the oxygen dis- 
placements: oxygens within a [011]c-oriented plane give 
raise to a dipole p along [Olljc. It is apparent that the 
addition of all such dipoles gives no net polarization in 
this case, a result that can be viewed as a consequence 
of the three-dimensional nature of the Og network. 

In conclusion, we have found that the strongest struc- 
tural instability in LTO's high-symmetry phase is a very 
common and simple one: it involves concerted rotations 
of the Og octahedra, much alike the AFD modes respon- 
sible for the structural phase transitions in most ABO3 
crystals with the ideal perovskite structure. AFD dis- 
tortions in the usual perovskites are well-known to be 
non-polar, a feature that can be viewed as a consequence 
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TABLE II: Energy difference between the FE and PE phases 
of La2Ti207 and BaTiOa. The FE phases are considered 
under several elastic constraints (see text), which we denote 
"fully-relaxed" (no constraint), "PE volume", and "PE cell". 
For both LTO and BTO, the energies are normalized to the 
volume of the PE phase so that a comparison between com- 
pounds can be made. Results given in meV/A^. 





LazTiaOy 


BaTiOs 


fully-relaxed 


-1.097 


-0.110 


PE volume 


-1.042 


-0.099 


PE cell 


-0.938 


-0.051 



of the symmetry and three-dimensional character of the 
lattice. In contrast, because the structure of LTO is split 
in layers comprising an even number of perovskite planes, 
the Og-rotation mode gives raise to a net polarization in 
this case. Hence, since the occurrence of a spontaneous 
polarization in LTO relies on the layered topology of the 
lattice, it seems appropriate to describe this compound 
as a topological ferroelectric. Let us stress that LTO's 
peculiar lattice topology does not seem essential for the 
structural instability to exist, but it is critical for it to 
have a polar character. 

C. Nature of the Cmcm-to-C'mc2i transition 

In addition to the already discussed topological nature 
of the leading FE instability, the transformation between 
LTO's Cmcm and Cmc2i phases presents a number of 
interesting aspects that we discuss in the following. Most 
importantly, our results show that the transition is driven 
by the combined action of the two soft modes discussed in 
the previous Section, and suggest that such a cooperation 
is critical for it to occur at a very high temperature. 

To study a structural phase transition quantitatively, 
one can start by comparing the energies of the phases 
involved. Table II shows our results for LTO, along with 
the analogous data for the Pmim (PE) and PAmm (FE) 
phases of BTO. Note that the energy change per unit 
volume involved in LTO's transition is about one or- 
der of magnitude greater than the corresponding one for 
BTO. Such an enormous difference is compatible with 
the experimentally measured Curie temperatures, which 
are 1770 K and 400 K for LTO and BTO, respectively. 
Hence, our first-principles results for the energetics of 
LTO's FE transition seem consistent with experiments. 

Table II also contains information about the elastic 
deformation that accompanies the FE transitions. The 
properties of FE oxides with the ideal perovskite struc- 
ture are known to be strongly sensitive to cell strains. 
This is clearly reflected in the results for BTO: If the 
compound is forced to keep the cubic cell of the PE struc- 
ture, the energy gain involved in the FE transition is 
reduced by half (i.e., it drops from —0.110 meV/A'^ to 
—0.051 meV/A^). On the other hand, if we only im- 



pose the PE volume (so that the FE cell is allowed to 
deform and acquire a. c/a ^ 1 ratio), the energetics of 
the transformation is not strongly affected. In contrast, 
the analogous elastic constraints have a relatively small 
effect in the case of LTO (i.e., the energy gain drops from 
-1.097 meV/A^ to -0.938 meV/A^ when we impose the 
PE cell). Such a weak coupling between strain and the 
FE distortion is probably reflecting the AFD-like char- 
acter of the instability)^ this result also suggests that 
LTO will display relatively small piezoelectric effects as 
compared with regular FE perovskites. 

Let us now consider the atomic displacements that 
characterize the Cmcm-to-Cmc2i transformation. In 
most materials undergoing displacive transitions, the 
atomic distortion connecting the HS and LS phases is 
essentially captured by the eigenvector of the instabil- 
ity mode that triggers the transformation. One can eas- 
ily quantify this by constructing a distortion vector AX 
comprising all the atomic displacements associated to the 
HS-to-LS transition, and expressing it in the basis formed 
by the eigenvectors of the force-constant matrix of the 
HS phase. (For simplicity, we constructed our distortion 
vectors AX using the atomic positions of a LS phase that 
is forced to have the same unit cell as the HS phase.) For 
BTO, this analysis led us to the expected result: the soft 
mode of the PE phase captures 98% of the atomic distor- 
tion involved in the PmSm-to-Pimm transition. How- 
ever, the result for LTO was qualitatively different: the 
strongest instability mode (^j, with ki = —0.81 eV/A^) 
captures 81% of the total distortion, and the second in- 
stability mode (^21 with K2 = —0.01 eV/A^) contributes 
with a 15%. (For both compounds, no other mode con- 
tributes more than a 1%.) Hence, our results suggest that 
the two soft modes and ^2 pl^^Y ^tn important role in 
LTO's FE phase transition. 

Such a large contribution of ^2 to the structural trans- 
formation seems incompatible with the computed mode 
stiffnesses; indeed, the values of ki and K2 would suggest 
that ^2 is about 80 times weaker than as an instabil- 
ity. To resolve this apparent contradiction, we computed 
how LTO's energy changes when the Cmcm phase is dis- 
torted according to the individual and ^2 modes. As 
shown in Fig. 4(a), we found a large energy reduction 
associated with the distortion, while the ^2 instabil- 
ity is almost negligible. We also computed the energy 
variation as the total distortion AX is frozen in. Re- 
markably, as compared with the results for ^j^, the AX 
curve in Fig. 4(a) presents a much deeper minimum cor- 
responding to a much larger distortion amplitude. This 
is a new indication that the soft mode cannot explain 
LTO's FE transition by itself. This result contrasts with 
the situation for BTO [Fig. 4(b)], where the FE soft mode 
captures the energetics of the structural transformation 
almost exactly. 

The results of Fig. 4(a) suggest that there is a strong 
and cooperative coupling between the two instability 
modes of LTO. To better describe this effect, let us write 
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FIG. 4: (Color on-line.) Panel (a): Variation of the energy of 
La2Ti207 as the Cmcm phase is distorted according to several 
atomic displacement patterns, namely, those corresponding to 
the soft modes and well as the distortion AX that 

connects the Cmcm and Cmc2\ phases. The 44-atom cell of 
the Cmcm phase is kept fixed in all calculations. Panel (b): 
Same as panel (a) for the PE {Pm3m) to FE {P4mm) tran- 
sition of BaTiOs. The inset shows the atomic displacements 
associated to the single FE soft mode of BTO (Ba atoms at 
the corners of the cubic cell; Ti atom at the center of the Oe 
octahedron). For both LTO and BTO, the energies are nor- 
malized to the volume of the PE phase so that a comparison 
between compounds can be made. 



the energy of the crystal as the following Taylor series: 

I / 3 I II 2 2 , III 3 I /-o/ 6\ 

where jg energy of the Cmcm phase, and ui 

and U2 are, respectively, the amplitudes (in Angstrom) of 
the $,1 and ^2 distortions. Note that this expression for 
the energy is greatly simplified by symmetry, and that 
the quadratic parameters coincide with the mode stiff- 
nesses. We have included in Eq. ^ only the lowest- 
order couplings between ui and U2, which are quantified 
by the primed 7 parameters. By fitting to the $,1 and 
^2 energy curves of Fig. 4(a), we got ki = —0.79 eV/A^ 
and K2 — 0.00 eV/A^, in fair agreement with the val- 
ues obtained from the diagonalization of K] we also got 
ai = 0.65 eV/A" and as = 0.66 eV/A^. Then, to fit 
the AX curve we considered distortions characterized 
by Ui/u2 = 0.81/0.15, as it corresponds to the Cmc2i 
phase. Note that, in this case, the fitted quartic coef- 
ficient is a combination of the a and 7 parameters of 
Eq. Given the large ui/u2 ratio associated with the 
AX distortion, and in view of further tests discussed in 
Section III.D, it seems reasonable to assume that 7' dom- 



inates over 7" and 7"'; we thus got 7' = —0.35 eV/A*. As 
compared with the computed a coefficents, this clearly is 
a very strong anharmonic couphng that favors the com- 
bined + ^2 distortion!^ 

To the best of our knowledge, such a cooperation be- 
tween soft modes is rare among perovskite oxides. There 
are many examples of materials in which several strong 
instabilities exist; in most cases, the strongest one leads 
to a phase transition that tends to suppress, partially or 
totally, the other instabilities. The competition between 
the FE and AFD soft modes in compounds like SrTiOs is 
a representative and well studied caseM- In contrast, we 
find that the reciprocal enhancement of the two soft FE 
modes is critical to explain the structural transformation 
in LTO. The mode is clearly the leading instability, 
and it would occur even in absence of ^2- Yet, as the re- 
sults in Fig. 4(a) show, the magnitude and strength of the 
transformation are boosted by the interaction between 
and ^2- Thus, our results seem to suggest that LTO owes 
its very high Tc to such an interaction; indeed, in view of 
Fig. 4(a), it is hard to imagine LTO's Tc would remain 
essentially the same if the coupling between and ^2 
was suppressed. 

Studying from first-principles the temperature depen- 
dence of the and ^2 distortions, and the details of how 
these two soft modes freeze in at LTO's FE transition, 
is a challenging endeavor that remains for future work. 
Nevertheless, a few observations can be made based on 
the present results. In principle, one could try to ap- 
proach the problem by introducing a Landau potential 
of the form: 



F-Fo= ^A,{T - Tc)Ql + \BiQi- 



1 

2' 

C'QlQ2 



\b2Q' 



2~ 

2/o2 



(3) 



C'QlQl + C 



'Q1Q2 : 



which is written in terms of two one-dimensional order 
parameters, Qi and Q2, that have the same symmetry. In 
the following heuristic argument we will consider only the 
Q1Q2 crossed term - in accordance with our above con- 
jecture regarding the couplings in Eq. and because 
this is the most relevant crossed term in vicinity of the 
phase transition-, ■^^ thus assuming that C" = C" = 0. 

If we were dealing with a simple transition, we would 
have a primary order parameter Qi corresponding to 
the unstable eigenmode of the K matrix of the high- 
symmetry phase; the temperature dependence of the 
Landau potential would be restricted to the term, 
as indicated in Eq. and we would have positive Ai 
and Bi coefficients. Further, a secondary order param- 
eter Q2 would be stable by itself, with positive A2 and 
B2 coefficients. Hence, the corresponding Landau theory 
would predict a phase transition at T — Tc, with 



1/2 



(4) 
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below the transition temperature. (The indicated formu- 
las were derived assuming that we remain close to the 
transition temperature.) Then, Q2 would present the 
following temperature dependence below Tc: 



01 



3/2 



(5) 



In a simple case, the coefficients in Eq. ([2]) may be a good 
approximation to the coefficients in Eq. ([3]) in the limit of 
low temperatures. Thus, by supplementing the Eq. ^ 
deduced from first-principles with a piece experimental 
information - i.e., the value of the transition tempera- 
ture Tc -, we could construct the corresponding Lan- 
dau potential and obtain a quantitative description of 
the temperature dependence of Qi and Q2- 

The doubts quickly appear when one tries to apply 
this model to LTO. It seems natural to associate Qi and 
Q2 with our computed FE soft modes and respec- 
tively. However, that identification implies that the tran- 
sition temperature for Qi is independent of Q2] such an 
approximation would be an awkward one, given the large 
influence that ^2 has in the energetics of the Cmcm-io- 
Cmc2i transformation. Additionally, we would be forced 
to speculate regarding the T-dependence of the A2 coef- 
ficient, as our guess for this parameter - i.e., the K2 of 
Eq. ([2|) which was found to be essentially zero - does not 
comply with the usual requirements for the quadratic 
coefficient of a secondary mode. Note that, for exam- 
ple, if we had a small value of A2 and thus a dominant 
B2 term, the behavior with temperature of both Qi and 
Q2 would vary: Qi would have the same functional T- 
dependence but with a different prefactor, and Q2 would 
go as {Tc - T)i/2 instead of (Tc - r)3/2. The resuhs 
of this study do not allow us to resolve such details of 
LTO's transition. 



D. Polarization and response properties 

Tabic III shows the computed spontaneous polariza- 
tion, lattice-mediated dielectric tensor, and piezoelectric 
tensor for the Cmc2i phase of LTO. For comparison, we 
also show results from the literature corresponding to the 
P2i phase that is stable at room temperature. 

We obtained a value of 0.29 C/m^ for the sponta- 
neous polarization, which is comparable with the result 
of 0.38 C/m^ that we obtained for the tetragonal phase of 
BTO. The dielectric response is also very significant, as 
the obtained values are comparable with those typical of 
FE perovskites (e.g., we got ezz = 23 for the z-polarized 
tetragonal phase of BTO). As regards piezoelectricity, the 
computed responses are considerable but not particularly 
large; for example, for the low-temperature rhombohe- 
dral phase of BTO, the piezoelectric coefficients reach 
values of 200 pC/N)^ while for LTO we got maximum 
values of about 20 pC/N. LTO's relatively small piezo- 
electric response seems compatible with the minor role 



TABLE IIL Non-zero components of the spontaneous polar- 
ization {P^ , given in C/m^), lattice-mediated dielectric ten- 
sor (ea/s), and piezoelectric tensor (dai, where I labels strain 
components in Voig notation, given in pC/N) of the Cmc2\ 
phase of La2Ti207. For the dielectric tensor, the clamped- 
cell response is given in parenthesis.— For the piezoelectric 
tensor, the clamped-ion response is given in parenthesis. We 
also show experimentalii and theoretical^^ values for the P2i 
phase of LTO that is stable at room temperature (Troom). 
Note that the results in Refs. [3 and [l^ have been adapted 
to our choice of Cartesian a:xes. 





Cmc2\ phase 


P2i phase 


this work 


exp. (Troom)^ theoryi^ 




0.29 
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that the cell strains play in determining the energetics of 
the CmcTO-to-Cmc2i phase transition, as mentioned in 
the discussion of Table II. 

The large spontaneous polarization obtained may seem 
incompatible with the AFD-like character of the struc- 
tural instability {^i) that dominates the Cmcm-to- 
Cmc2i transition. To clarify this point, we performed al- 
ternative calculations using linear-response expressions^ 
based on the Born effective-charge tensors Z* that quan- 
tify the polarization change associated to the displace- 
ment of an individual atom i. 

The polarization reported in Table III was obtained 
in the standard way: We used the Berry-phase theory of 
King-Smith and Vanderbilt^ to compute the variation of 
P as the Cmc2i structure is deformed into a symmetry- 
equivalent one with opposite polar distortion; then, 
is half of the computed polarization change. P^ can also 
be estimated using the approximate formula 



pti 



i/3 



AX 



1/3 



(6) 



where AX is the vector capturing the distortion that 
connects the Cmcm and Ctoc2i phases, i labels the 
atoms in the unit cell, and a and /3 label spatial direc- 
tions. Using different choices for the effective-charge ten- 
sors (i.e., those computed for the Cmcm phase, as well 
as the corresponding results for the Cmc2i phase subject 
to different cell constraints) we obtained values of in 
the 0.25-0.32 C/m^ range, which are perfectly compatible 
with the result in Table III. 

Then, we used the effective-charge tensors of the 
Cmcm phase (given in Table IV) to compute the polar- 
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TABLE IV: Computed effective-charge tensors (given in units 
of elementary charge) for the Cmcm (PE) and Cmc2\ (FE) 
phases of La2Ti207. To facihtate the comparison between 
phases, we list the tensors corresponding to all the symmetry- 
inequivalent atoms of the FE phase. Atoms are labeled as in 
Table I. 



La(l) 



La(2) 



Ti(l) 



Ti(2) 



0(1) 



0(2) 



0(3) 



0(4) 



0(5) 



0(6) 



0(7) 



Cmcm (PE) 



Cmc2i 




ization change associated with the condensation of the 
and ^2 soft modes that dominate the Cmcm-io-Cmc2i 
transformation. The polar character of the modes is 
quantified in terms of the mode effective charges: 



If) 



(7) 



TABLE V: Computed effective-charge tensors (given in units 
of elementary charge) for the Pm3m (PE) and Pimm (FE) 
phases of BaTiOs. The tensors are given in the conven- 
tional Cartesian axes for a rectangular lattice. They corre- 
spond to the following atoms (PE phase positions given in 
relative units): Ba located at (0,0,0), Ti at (1/2,1/2,1/2), 
0(1) at (1/2,1/2,0), and 0(2) at (0,1/2,1/2). 0(1) and 0(2) 
are symmetry-equivalent in the PE phase, but not in the z- 
polarized FE phase. 



PmSm (PE) 



Pimm (FE) 




We obtained Zi^z — l-8e and Z2^z = 12. Oe, where e is the 
elemental charge. These results confirm that the is 
weakly polar, in accordance with its AFD-like nature. In 



contrast, the second soft mode ^2 is found to have a con- 
siderably polar character. It can then be trivially shown 
that the two soft modes have a very similar contribution 
to the total spontaneous polarization of the Cmc2i phase 
of LTO, in spite of the fact that embodies 81% of the 
PE-to-FE distortion. Hence, the two-mode character of 
the transition is critical to explain the relatively large 
of the Cmc2i phase. 

Unfortunately, we were unable to find experimental re- 
sults for the ferroelectric and response properties of the 
Cmc2i phase of LTO. Table III shows some results for 
the P2i phase of the compound that is stable at room 
temperature. Interestingly, the dielectric and piezoelec- 
tric responses measured for this phase seem compatible 
(at least in magnitude) with our values for Ctoc2i. As 
regards the spontaneous polarization, the value for the 
P2i phase is smaller than ours by a factor of 6; this ex- 
perimental result was essentially confirmed by the first- 
principles study of Ref. [l^, which employed a DFT the- 
ory similar to ours.^ Hence, we can tentatively conclude 
that the Cmc2i-to-P2i transformation, which is experi- 
mentally determined to occur at 1053 involves a re- 
duction of LTO's spontaneous polarization. While such 
a reduction in the magnitude of as the temperature 
decreases is not typical in FE perovskites, in principle 
there is no reason to question this possibility. 

Finally, let us discuss an intriguing possibility sug- 
gested by the very different magnitudes of the computed 
mode charges Zi^z and Z2.z- Since quantifies the cou- 
pling between a polar mode and an applied electric field, 
we can expect ^2 to be much more reactive to an external 
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bias than can thus imagine the following pos- 

sibility: to apply an electric field to LTO in its Cmc2i 
phase and switch only the pari of the spontaneous polar- 
ization associated to More precisely, if {u\,u\) repre- 
sents the FE phase discussed so far, an electric field might 
allow us to take the material into a qualitatively differ- 
ent polarization state {uY,u^) ~ {u\, —u\). LTO would 
thus be a four-state ferroelectric, as the {—u\, —u\) and 
{—uY, —u^2) variants would be accessible as well. A nec- 
essary condition for such a partial switching to occur is 
that the {u^, u^) state be a minimum of the energy. More 
precisely, we need the coupling between ^-^ and ^2 to be 
dominated by the 7"u^U2 term of Eq. ([2]): Note that, 
for \'y"\ much greater than |7'| and the four states 

{ztu\,:Lu\) would have essentially the same energy; in 
contrast, if 7' or 7"' dominate, the "^2"Switched" state 
would have a much higher energy and might not be a 
minimum of E[ui,U2)- In order to confirm or disprove 
such a partial switching, we considered the equilibrium 
structure of the Cmc2i phase (i.e., the one reported in 
Table I) and generated configurations in which the ^2 
distortion was inverted and given various magnitudes; we 
relaxed such transformed structures and invariantly ob- 
tained the original Cmc2i phase as final result. Hence, 
while we did not explore the E(ui,U2) energy landscape 
in detail, our results clearly indicate there is no stable 
"^2-switched" state. Equivalently, this implies that the 
7' and 7"' terms of Eq. ([2]) dominate over 7" (which sup- 
ports the assumption made in Section III.C). 



E. BaTiOs-like ferroelectric modes in La2Ti207 

Thus far we have discussed the main features of LTO's 
high-temperature FE transition, showing that this com- 
pound is very different from the FE oxides with the ideal 
perovskite structure. In particular, our results seem to 
suggest that LTO and BTO have little in common, in 
spite of the fact that they share the same building block: 
TiOe octahedra with Ti''+ in the electronic configu- 
ration. In the following we show that such a conclusion 
would be a deceptive one. 

Ferroelectricity in the usual FE perovskite oxides re- 
lies on strong dipole-dipole interactions whose fingerprint 
is the anomalously large magnitude of the Born effec- 
tive charges of some atomic species4^ A prototypical ex- 
ample of this behavior is BTO, for which we computed 
the effective-charge tensors shown in Table V. Most no- 
tably, in the PE phase the Ti^+ cations display Z* val- 
ues above 7e, almost doubling their nominal ionization 
charge. Analogously, the displacement of the O atoms 
towards the Ti has a dynamical charge of — 5.77e asso- 
ciated to it, almost tripling the nominal value of — 2e. 
It is well-known that this effect is related with a partial 
covalent character of the Ti-0 bond, and an increased 
hybridization of the 0-2p and Ti-3d orbitals as the Ti 
and O atoms approach. Thus, given this chemical ori- 
gin, the effective charges become less anomalous once the 



FE distortion freezes in: according to the results in Ta- 
ble V, Ti's charge of 7.31e gets reduced to 5.83e, while 
O's charge of -5.77e falls to -4.79e. 

As apparent from Table IV, some Ti and O atoms 
in LTO also present anomalously large effective charges 
that reach values similar to those obtained for BTO. In 
the case of the Ti atoms, the computed Z* tensors are 
rather isotropic, with the maximum Z* values [i.e., 6.91e 
for Ti(l) and 7.45e for Ti(2)] corresponding to displace- 
ments along the in-layer directions x and z. In the case 
of the O atoms, 0(6) and 0(7) clearly present the largest 
values, which exceed — 5e as in BTO. Such giant dynam- 
ical charges corresponds to displacements along the x di- 
rection, for which LTO presents infinite chains of TiOe 
octahedra (see Figs. 1 and 2) as those in the ideal per- 
ovskite structure. Hence, as regards their polarizability 
properties, the results for the Ti and O atoms in LTO are 
strongly reminiscent of the behavior that is well-known 
for BTO. 

Do such anomalous effective charges lead to FE insta- 
bilities in LTO? In accordance with the above discus- 
sion, the answer to this question is a negative one: the 
Cmcm phase presents only two F-point instabilities (i.e., 
the already discussed and ^2) the Cmc2i is sta- 
ble against F-point perturbations^^. Nevertheless, by in- 
specting the if-eigenmodes of the Cmcm phase, it is easy 
to find several low-energy FE modes that are BTO-like, 
i.e., they involve the displacement of the Ti atoms away 
from the center of nearly undistorted Og octahedra4^i^ 
For example, we found a marginally stable and strongly 
polar mode for which we computed Kg = 0.26 eV/A^ and 
Zs^x = 21. 7e (this distortion is x-polarized and has B^u 
symmetry) ; this mode gives an enormous contribution to 
the exx dielectric response:>iS, the obtained value exceeds 
600. Interestingly, these BTO-like modes become stiffer 
in the Cmc2i phase; thus, they do not give raise to any 
anomalously large contribution to the dielectric response 
reported in Table III. 

In conclusion, our analysis shows that LTO presents 
obvious traces of the FE instabilities of BTO. Such a 
transferability of instabilities was demonstrated by one of 
us in an hexagonal polymorph of BTO. "^^ In that case, the 
FE phase transition was shown to be driven by soft modes 
that are an almost perfect match of the usual BTO-like 
FE instability. In LTO, such modes are very low in en- 
ergy, but still stable; further, they clearly compete with 
the dominating AFD-like instability, and become stiffer 
once the Cmcm-to-Cmc2i transition occurs. Neverthe- 
less, noting that BTO-like FE distortions tend to be very 
sensitive to cell strains, our results suggest the possibil- 
ity that such instabilities might be induced in LTO by 
suitable strain engineering or chemical substitution, or 
that they might occur spontaneously in similar layered 
perovskites. 
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F. Implications for work on magnetoelectrics 

Finding materials that display large magnetoelectric 
(ME) effects (i.e., a large magnetic reaction to an ap- 
plied electric field) at room temperature is a major chal- 
lenge that remains to be successfully tackled. The dif- 
ficulties involved in the design of good magnetoelectrics 
have been discussed elsewhere At present, the strate- 
gies that seem most promising rely on finding systems 
that satisfy the following two conditions: (1) their atomic 
structure must react strongly to an applied electric field 
(i.e., we are looking for good dielectrics), and (2) the 
field-induced distortions must have a large effect on the 
magnetic interactions (as emphasized in Refs. 48 andlioh. 

While there are well-known strategies to comply with 
the first condition,''^ satisfying the second one is prov- 
ing much more difficult. In fact, the existing quantita- 
tive studies indicate that the dielectric response of ME 
multiferroics like BiFeOa is dominated by modes that 
have small magnetostructural couplings associated to 
them. 5- It is thus interesting to consider the alternative 
approach adopted by Benedek and Fennie>^ These au- 
thors noted that modes involving Oq rotations are likely 
to be strongly coupled with the magnetism of perovskite 
oxides, as such AFD-like distortions usually control the 
nature and magnitude of the main magnetic interac- 
tions (e.g., the metal-oxygen-metal super-exchange and 
Dzyaloshinskii-Moriya couplings). Hence, they looked 
for materials in which AFD-like distortions are related 
with (or lead to) ferroelectricity, as in such cases one 
could use an electric field to act upon the Oq rotations. 
Note that the seeked connection between AFD-like dis- 
tortions and ferroelectricity is an exotic one, since the 
Og-rotational modes are strictly non-polar in compounds 
with the ideal perovskite structure (see Section III.B). 
Accordingly, ferroelectricity-related Og-rotational modes 
have been found in non-ideal perovskites, e.g., in the 
PbTiOa/SrTiOs artificial superlattices^ and the layered 
compound Ca3Mn2 07.^'^ In both cases, ferroelectricity 
has an improper character, and a combination of various 
AFD-like modes must occur for an spontaneous polariza- 
tion to appear. 

There is no need to emphasize the importance of our 
findings for LTO in the context of magnetoelectrics. In 
LTO's case, the FE soft mode is AFD-like already; hence, 
it can freeze in and give raise to a sizeable spontaneous 
polarization by itself, without the need of any accom- 
panying distortion. Further, LTO's Og rotations couple 
directly (bi-linearly) with an applied electric field, which 
might prove advantageous for the purpose of obtaining 
large ME effects. (In improper ferroelectrics, the cou- 
pling between an applied field and the AFD-like modes 
will typically be a higher-order effect.) Indeed, the proper 
ferroelectricity driven by Oq rotations that occurs in LTO 
seems to be the dreamed FE instability from the view- 
point of ME applications. 

One would thus like to obtain LTO-like ferroelectric- 
ity in a magnetic material. The most obvious possibil- 



ity is to consider the substitution of Ti by Mn in LTO, 
so as to form La2Mn207, a crystal that we have not 
found described in the literature. In La2Mn2 07 we would 
have manganese in the Mn^+ ionization state, most likely 
in the high spin ifg^g electronic configuration; thus, we 
can expect this crystal to be insulating. Further, since 
the Ti-0 chemistry does not seem to play any relevant 
role in LTO's FE instability, we can expect La2Mn207 
to present a FE transition analogous to LTO's. Hence, 
this material seems an excellent candidate to satisfy the 
condition (2) mentioned above and thus display large 
ME effects. Additionally, one would like the FE tran- 
sition to occur near room temperature, so as to benefit 
from the enhancement of the functional responses near 
Tc [in the spirit of the condition (1) mentioned above]. 
In this sense, it may be useful to note that materials like 
Sr2Ta207^ or Sr2Nb207^ present FE transitions that 
seem similar to LTO's but occur at lower Tc's: 166 K 
and 1615 K, respectively. This suggests that exploring 
alternative compositions is a promising route to tunc the 
temperature of the Oe-rotational FE transition. 



IV. SUMMARY AND CONCLUSIONS 

We used first-principles methods to study the origin 
of ferroelectricity in the layered perovskite La2Ti207 
(LTO), one of the materials with highest Curie temper- 
ature known (Tc — 1770 K). To do so, we carried out 
for LTO a research program that has been repeatedly 
and successfully applied to the investigation of displacive 
phase transitions in ferroelectric (FE) oxides with the 
perovskite structure. Our results allowed us to charac- 
terize LTO's high-temperature FE transition, which was 
found to present very peculiar features. 

We found that ferroelectricity in LTO is very differ- 
ent from what occurs in the well-known FE oxides with 
the ideal perovskite structure, such as BaTiOs (BTO) or 
PbZri_:rTij;03 (PZT). Indeed, the dominant FE instabil- 
ity of this compound has little in common with the text- 
book picture of positive charges moving against negative 
charges in an ionic insulator; instead, it involves con- 
certed rotations of the oxygen octahedra that form the 
perovskite framework. Hence, LTO's high-temperature 
structural transition is reminiscent of the behavior of 
simple perovskite oxides, but not the FE ones: LTO's FE 
distortion is much alike the Og rotational soft modes that 
drive the anti-ferrodistortive phase transitions of SrTiOa, 
LaAlOs, LaMnOs, and many other non-polar perovskite 
crystals. 

Hence, we found that the existence of ferroelectricity in 
LTO at record-high temperatures is the result of struc- 
tural distortions as those occurring in many perovskite 
oxides that are not ferroelectric. As discussed here, the 
solution to this puzzle has to do with LTO's layered struc- 
ture: Because the lattice of oxygen octahedra is trun- 
cated in this compound, the Og-rotations acquire a polar 
character and give raise to a macroscopic polarization. 
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We can thus describe LTO as a topological ferroelectric, 
since it owes its spontaneous polarization to the layered 
topology of its structure. Note that such a peculiar topol- 
ogy does not seem essential for LTO's high-temperature 
transition to occur, but it is critical for it to have a fer- 
roelectric character. 

We quantified energetics of LTO's FE transformation, 
the results being consistent with the very high temper- 
ature at which it happens. Most remarkably, we found 
that the structural distortion connecting the PE (Cmcm) 
and FE (Cmc2i) phases presents large contributions 
from two modes, namely, the above-mentioned strong in- 
stability consisting of Oq rotations, and an isosymmet- 
ric weakly-unstable mode that involves deformations of 
the oxygen octahedra. Further, we found that the large 
energy change associated to the C'mcm-to-Cmc2i tran- 
sition depends crucially on the simultaneous occurrence 
of both modes. We were able to describe such effects in 
terms of a very strong and cooperative anharmonic cou- 
pling, and briefly discussed the possibility of constructing 
a phenomenological theory of such a two-mode transition. 
Interestingly, LTO's behavior is in strong contrast with 
what (to the best of our knowledge) is most common 
among perovskite oxides, where structural instabilities 
tend to compete and suppress each other. 

We calculated the polarization and response proper- 
ties of LTO's FE phase, obtaining results that are con- 
sistent with related experimental information. Interest- 
ingly, the computed electric properties revealed clear sim- 
ilarities between LTO and prototype ferroelectric BTO. 
For example, we obtained anomalously large Born effec- 
tive charges for some Ti and O atoms in LTO, a fea- 
ture that is known to be the fingerprint of the FE in- 
stabilities in compounds like BTO and PZT. Further, 
our results showed that LTO presents nearly-unstable 
modes that are strongly polar and involve atomic distor- 



tions that essentially match BTO's FE instabilities. Such 
traces of BTO-like behavior in LTO suggest the intrigu- 
ing possibility that conventional ferroelectricity might be 
induced in this compound upon suitable modifications 
(e.g., strain engineering or chemical substitutions), or 
that such a behavior might occur spontaneously in simi- 
lar materials. 

Finally, let us emphasize the implications that our 
findings have for the design of new magnetoelectric 
multiferroics. In the context of magnetoelectric per- 
ovskite oxides, it would be ideal to have proper ferro- 
electricity driven by Og-rotational modes, so that an 
electric field can be used to tune the structural dis- 
tortions (e.g., metal-oxygen-metal angles) that control 
the main magnetic interactions (i.e., super-exchange and 
Dzyaloshinskii-Moriya). That is exactly the kind of 
ferroelectricity that we have found in LTO. Here we 
have briefly discussed how to obtain such an effect in 
a magnetic perovskite, proposing La2Mn2 07 as the most 
promising candidtate material. 

In conclusion, our theoretical study of La2Ti207 has 
revealed a wealth of interesting effects, some of which 
may have important implications for current work on 
multifunctional oxides. We thus hope this study will 
stimulate further investigation of these layered per- 
ovskites and the novel functionalities that they may offer. 
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